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incompressible, viscous, and fully developed. The nonlinear 
equations under suitable boundary conditions were scaled to 
a system of dimensionless PDE, which were reduced further 
to ordinary differential equations using the perturbation 
conditions. The perturbed system of ordinary differential 
equations was solved using the Laplace method. The blood 
flow profiles, such as velocity, volumetric flow rate, lipid 
concentration profiles such as lipid concentration and rate of 
mass of the lipid transfer, and temperature profiles such as 
blood temperature and the rate of heat transfer are obtained 
and the effects of the pertinent parameters such as magnetic 
field, radiation absorption, and Brinkman Number were 
discussed. 
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The results show that the velocity increases with an increase in Brinkman number, Grashof number, solutal 
Grashof number and the porosity, while it decreases with an increase in Schmidt number, Prandtl number, 
radiation absorption, chemical reaction, and magnetic field. From the investigation, we could deduce that 
there is an increase in blood temperature with the increase in Brinkman number and a decrease in 
temperature due to an increase in Schmidt number, Prandtl number, radiation absorption, chemical reaction, 
and oscillatory frequency. In conclusion, we have been able to formulate a mathematical representation of 
blood flow through a sclerotic microchannel, obtained an exact solution to the problem, and presented 
results that might be useful for mathematicians and clinicians. 


1 Introduction 


Blood, the heart, and blood vessels are all part of the cardiovascular system. Blood is essential because it 
serves as a transport agent in the human body. Unfortunately, plaques in human blood vessels such as 
arteries and capillaries can disrupt normal blood flow, leading to cardiovascular diseases such as heart attack 
and stroke. Because of the implications for medicine and fluid mechanics, abnormal blood flow has piqued 
the interest of many researchers. Blood can be classified as a non-Newtonian fluid, so studies that involve 
modeling blood flow should take this into account. Humans are constantly confronted with challenges that 
disrupt the normal circulation of blood. This includes, among other things, physical exercises, vehicle travel, 
and application, according to Bunonyo et al. [1]. 


Over the years, lots of mathematicians and scientists have done ample research on the flow of blood in the 
human circulatory system and have discovered some useful results; here are a few of those studies 
supporting this research work. Bunonyo et al. [2] studied the impact of treatment parameters on blood flow 
in an atherosclerotic artery. Bhatti et al. [3] investigated the heat transfer analysis of peristaltic induced 
motion of particle fluid suspension with variable viscosity: the clot blood model. The influence of blood 
flow in large vessels on temperature distribution in hyperthermia was analyzed by Lagendijk [4]. Srivastava 
[5] discussed the analysis of the flow characteristics of blood flowing through an inclined tapered porous 
artery with mild stenosis under the influence of an inclined magnetic field. Bunonyo and Amos [6] 
investigated the effects of treatment and radiation on oscillatory blood flow through a stenosed artery. The 
study of slip effects on the unsteady MHD pulsating blood flow through a porous medium in an artery under 
the effect of body acceleration was carried out by Eldesoky [7]. Makinde and Mhone [8] investigated the 
heat transfer to MHD oscillatory flow in a channel filled with porous medium. Unsteady heat transfer to 
oscillatory flow through a porous medium under a slip condition was explained by Hamza et al. [9]. 
Choudhury and Das [10] studied the heat transfer to MHD oscillatory viscoelastic flow in a channel filled 
with porous medium. Biswas and Chakraborty [11] gave the pulsatile flow of blood in a constricted artery 
with body acceleration. Hanvey et al. [12] developed a model of heat and mass transfer in the oscillatory 
flow of a non-Newtonian fluid between two inclined porous plates placed in a magnetic field. Plourde et al. 
[13] looked into alterations in blood flow through arteries following atherectomy and the impact on pressure 
variation and velocity. Mathematical modeling of blood flow through vertebral arteries with stenosis was 
studied by Ali et al. [14]. Mathur and Jain [15] discussed the mathematical model of non-Newtonian blood 
flow through an artery in the presence of stenosis. Ali and Asghar [16] studied the oscillation channel flow 
for non-Newtonian fluids. 


Varshney, Katiyar, and Kumar [17] investigated the effect of a magnetic field on the blood flow in an artery 
having multiple stenosis. Elangovan and Selvaraj [18] gave an insight into multiple stenosed arteries with 
periodic body acceleration in the presence of a magnetic field. Transport of MHD couple stress fluid through 
peristalsis in a porous medium under the influence of heat transfer and slip effects was investigated by 
Sankad and Nagathan [19]. Tripathi and Sharma [20] studied the effects of variable viscosity on MHD 
inclined arterial blood flow with chemical reaction. Modeling of arterial stenosis and its applications to 
blood flow was given by Pralhad and Schultz [21]. Bunonyo and Amadi [22] investigated the oscillatory 
flow of an electo-hydrodynamic fluid flow through a channel with a radiative heat and magnetic field. 
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Bunonyo and Eli [23] investigated the convective fluid flow through a sclerotic oscillatory artery in the 
presence of radiative heat and a magnetic field. 


In this research, radiation-brinkman number effects on blood flow in a porous atherosclerotic microchannel 
with the presence of a magnetic field, growth rate, and treatment are investigated using the Laplace method. 
It involves the use of mathematical models to represent blood momentum in an arterial channel, flowing 
through an atherosclerotic artery due to an exponential growth of cholesterol on the upper walls of the artery. 
The simulated results are presented graphically by investigating the varying pertinent parameters. 


2 Models Formulation 


We consider blood to be unsteady, electrically conducting, and incompressible viscous blood flowing 
through an atherosclerotic artery, which is presumed to be a cylindrical microchannel with a velocity 
w (rx), where r and x the directions of the flow. The flow in the azimuthal direction is considered to 
be zero. The flow is due to the pumping action of the heart; the pressure gradient is in the horizontal 


direction; and the magnetic field is perpendicularly applied to the direction of the flow of blood. The models 
governing the flow can be presented as: 


2.1 The geometry of atherosclerosis 


Ry [ eos a l at d Sx SÀ 
2 A 


R= (1) 


R, at O<x <d, 


where x = (a + z) (2) 


2.2 Blood Momentum Equation 


ôw OP" Ow 1dw) wo « T x : 
Pp or = ox ral je, aa 5 w -0oB;w + p,8B;(T -T,,)+ P,8Be(C =C.) 
(3) 
2.3 Energy Equation 
oT OT 1 ôT” i pas 
Poor ae = a F r r) Q, (T -T,)+Q; (c -C,) (4) 


2.4 Lipid Concentration Equation 


OC _ (25 „1 c) n(c -c.) 5) 


or” or? r or 
The corresponding boundary conditions are 


w =0,T =T 


w? 


C*=C, a r =R 
(6) 
w #0T +T, ,C* +C, ar =0 
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2.5 Dimensionless Scaling Parameters 


7 T, -T,,)Ry = : 
pera a w= WR, gr- ; )Ro de — eg oS 
A Ry 0 V, V, T-T V, 
p= Dkn í T, ae jem Q,Ry ,Rd, = Q,R; (C, -C,,) ,Gc = 8BAC, — Cn) Ry : 
vT, C, ~ Co Lp b kr (T, T) U, (7) 
2 at 3 _ 
M = BR, ae phi „S 7s „ô aau P Ps p A g=- Co 
pk k D, R, kn, A Mv, C=C, 
*2 2 
Q, = PW , Br = HW , ô 
U, (C, ~ C,) kr (T, -T,) Ry 
Applying the Dimensionless parameters in equation (7), equations (1)-(6) are reduced to 
R -rer (1+ cos2zx) at d,<x <1" (8) 
Pe li at OS x <d, 
1 A 
where x = —| d, +— 9 
” Al 0 5] a 
2 
OW Pa DOW | Tees tea Bea (10) 
Ot ôx \ er Fer k 
2 
poe a Ze ag Rd PrO + Bro (11) 
Ot ôr ror 
0d (p 10¢ 
Sc— =| —+-— |- Rd,Sc 12 
Ot E r or peo (1) 
The corresponding boundary conditions are 
w#0,040,¢40 at r=0 
(13) 


w=0,0=1,¢=1 at ie 
Ry 


3 Reduction to ODE by Perturbation Method 


We shall apply the perturbation method to reduce the partial differential equations (10)-(12) to ordinary 
differential equation, whence we can apply the appropriate method to solve the ODE. The solution to the 
PDE can be presented in the following form: 


iot 


w(r,t)=w,(r)e 
O(r,t)=6,(r)e 


g(r,t)=,(r)em (14) 
-2P _p gia 
Ox 


Using equation (14) in reducing the PDE, we have the ODEs as follows: 
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d’w, 1dw 

a? i i -piw = —F, —Gr@, -Gch (15) 
d°0, 1d0, 

ie Hor PaO = = -Bro (16) 
d’ ld 

h] h p h= (17) 


where 8? = 1 Mio „B: =(Rd, +ia) Pr, B; =(Rd,+io) Sc 
1 k 2 1 3 3 


The corresponding boundary conditions are 


w +0,0, #0,¢, #0 at r=0 
, , 18 
Wy = 0, 6, = e”, hy =e" a r= R ( ) 
R, 
4 Method of Solution 
The Laplace of the basic function can be stated as follows: 
L{w, (r)} =W (s)=[ my (re "dr (19a) 
0 
L{4,(r)}= 4, (s) =] (rJe "ar (19b) 
0 
L{h(r)}= (8) =] 4 (near (19c) 
0 
Restating equation (17), we have 
d'h 1 a 
Th Lah ip.) 4 =0 (194) 
We let 2}, = iĝ}, and then equation (19d) reduces to: 
Ph 1d | 
20 
Gt get Bath = (20) 
Applying Laplace method in solving equation (20), we state it as follows 
1 Shh {Sel ger {rg} =0 (21) 
r 


Simplifying equation (21), we have 


ndr Eh ladt) par {rh} =0=-2 (5% (5)-s (0) -4 (0))+ s4 (s)-8 „(0 J= Ba po 
a 


= (5% (s)— sp (0)-% (0)) +54, (s)—@ (0)- Ba e po = (23) 
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Ah, SA —=0 (24) 
ds (s + Ba) 
Simplifying equation (24), we have 


B 
h (s)= -= (25) 
(s + bzi ) 
Taking the inverse transform of equation (25), we have 
4 B, 
h (1) =L"§ ~——$—~ t = BJ, (Bur) = BJ (ibr) (26) 
(s + Ba) 
Note that J, (ir) =I, (A,r), so that equation (26) becomes 
(7) = Bly (Br) (27) 
Solving for the constant coefficient in equation (27) using the condition in equation (18), we have 
ei 
B =——_— (28) 
| 1,(Bh) 


Substituting equation (28) into equation (27), we have 


—iot 


$(r) = F] L (Br) (29) 


To obtain the lipid concentration profile, we shall substitute equation (29) into equation (14), which is 


g(r.) [aan (30) 


L (Bh) 
Retransforming equation (16), we have 
2 
cad a +A) 0, =—Brd, (31) 
dr’ E 


We let £,, =ip,, and then equation (31) reduces to: 


d’0, 1d0, 
Jr mH Pea 2 + pab = =-Brg, (32) 
Decompose equation (32) into homogenous and non-homogenous, which is 
dO, 1d0, 

+-+ + 86, =0 33 
dr? r dr + baa ( ) 
Solve equation (33) using Laplace method, which is 

d°0, dO 
ur “a +f Os + iL r6,}=0 (34) 


Simplifying equation (34), we have 
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d’6, db} n d, . 2 dO, 
ur 4 haf : |4-piL{r0,}=0=-2(s 0,(s)—30, (0) -4,(0)) +30,(s) -,(0)- 23 


r r 
(35) 

Simplifying equation (35), we have 
d0; S 

+ 0, (s)=0 (36) 

ds (s°+£;,) a(s) 
Solving equation (36), we have 
B 
O; (s) = > (37) 
(s + Pal 


Taking the inverse transform of equation (37), we have: 


5 B. E 1 ; 
O(a — 220: "sy = B,J, (B,,r) = B,J, (ibr) (38) 
[s + Ba) (s +a) 
Note that J, (i,r) =1,(,r), so that the homogenous solution of equation (32) is 


,(r) =BI, (Ar) (39) 


Recalling equation (16) after substituting the lipid concentration profile, we have 


a6, {1d0, m ei 
+—— — £50, =-Br| ——— |I r 40 

dr’ r dr B; 0 I, (Bh) o (b: ) ( ) 
Let the particular solution be in the form 
A, (r)=4 +4 (4,7) (41) 
Differentiate equation (41) twice according to the order of equation (40), which is 

dé, d*0, 
"a = A, BI, (Br), a = ABI (Br) (42) 
Substitute equations (42) and (41) into equation (40), we have 

1 o 
ABE (Br)+=(4 BI, (Br))- 2A - PA lo (Br)=-Br Iy (Br) (43) 
r I ( Bh) 
Simplifying equation (43) we have 
Br ei 

A, =0,A, =—| —— (44) 

b: í I 0 (Bh) ) 
Substitute equation (44) into equation (41), we have 

Br| e” 

O, (r)= =| — Mh (Ar (45) 

bol) Pa a) 


The general solution of equation (16) is the sum of equations (45) and (39), which is 
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0, (r) = Bl, (Br) + 5 Gh tA (46) 


Solving equation (46) using the boundary condition in equation (18), we have 
ei Brei@ 
B, =—— - > 
I (h) B:1,(Bh) 


Substitute equation (47) into equation (46), which is 


(47) 


—iøt 


a- rE ELED]! (Br) Al can) a a 


The temperature profile with the effect of the lipid concentration, we substitute equation (48) into equation 
(14), we have 


| i i Bla (bh)) ' fol ar)+ A L (bh) (8 ) : (49) 


The other key objectives of this research are to investigate impact of temperature and lipid concentration 
profiles on blood momentum. Let us recall the momentum equation as 


d wo, 1 dw, 
dr’ i dr 
Substitute equation (48) and (29) into equation (15), which is 


d°w 1 dw Gre” Br ei BrGr 
m r poh es fy fi a (Br J í Joes 2 Jalar) 


-piw = R,- Gr0, - Gch (50) 


Io (b-h) b: Io (Bh) B, 
(51) 
To solve the homogenous part of equation (51), we let it to be equal to zero, that is 
d “Wo 1 dw, 
= 52 
dr’ t dr Piw (52) 
We can further transform equation (52) as 
d'w, 1 dw, 2 
+—-— +(i w, =0 53 
dr? r dr ( iB, ) 0 ( ) 
We let f,, =if,, and then equation (53) reduces to: 
d’w. ldw. p 
+-—2 + Bw, =0 54 
d? rdr a e 
Solve equation (54) using Laplace method, which is 
d’w, dw A 
en s 


Simplifying equation (55), we have 
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rdr eih {Hah girir} =0=-2 (sws) -a (0)-8(0)) 4095 (8) (0) Z 


r ds 
(56) 

Simplifying equation (56), we have 
dw, s 

+ w,(s)=0 (57) 

ds (s? +p) ols) 
Solving equation (57), we have 
B 
Wo (s) = == (58) 
(s + B) 


Taking the inverse transform of equation (58), we have: 


w,(r)=L" pe = BL" TS = B,J, (Br) = B,J, (ifr) (59) 
Vs? +4) (5° + BA) 


Note that J, (ifr) =1, (A,r). so that equation (59) which is 
wy (1) = Bl (Br) (60) 
Recalling equation (51) after substituting the lipid concentration profile, we have 


d'w ldw p _ Gre” Br A ao BrGr r 
a rar a, a joe) gle p: Jen 


(61) 


The particular solution of equation (62) takes the form 

Wop (r)=A,+ Aly (Ar) + Alo (Br) (62) 

Differentiate equation (62) twice according to the order of equation (61), which is 

d’w, 3 
dr’ 


d 
Tor = (Ap (Byr) + ABL (Byr)) > 2 = (ABE ( Bar) + ABT (Br) (63) 


dr 


Substitute equations (63) and (62) into equation (61) and solving, we have 


P 1 | Gre Br 1 BrGr ae 
A, =47 =- l-7 | p4 = Ga] C 2 i 
; As B e b: } As ral a B, sz i 


Substitute equation (64) into equation (62), we have the particular solution as 
P 1 | Gre” Br 1 BrGr ee 
Wo, (r) =o + l1- l I, (Ar te l I (Br 
"xi ale E (Bh) AJ T Ha By Na ) 


The general solution of equation (62) is the sum of equations (65) and (60), which is 


(65) 
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P 1[( Gre Br 1 BrGr e`” 
w (r )= BI, (Br )+ B e- B 2) aes )+ p [ces B aoe 


(66) 
- R 1 (Gre (,_ Br E AE 
w(ra)=[ But oe a am f ay nee xe B (nt j 
(67) 
where B, = a (z | 5 R — = [ces a) 
P L(Bh) B P I, (Bh) BI, (Bh) B; 


The flow cab be calculated mathematically as 


r=h 


iot P, 1[ Gre“ _ Br al BrGr iar 
a Jeen Bi wi (Bh) |: B: Jets op [aes A lga] 
Simplifying equation (68), we have 


r=h r=h 1 Gre” Br 

B,rl,(f,r)dr+ 1 rI,(B,r)dr 

vase um” [eS (Soe 
= (d 


=] i ze) a | "Pr 
+ Gc+ rI (Br)dr+ | dr 
r=0 VA b: I) (Bh) i ( i ) i B 


After integration and substituting the limit of integration, we have 
Gre™ Br 
Bh rt ( J] Byh)+ 
a Al )+ B, Be 0 I, (Bh) b: ( : ) (70) 


h BrGr ei Ph 
Ge 1, (Bh 
Fel EB 5) (Ail) + 


The rate of lipid transfer can be calculated mathematically as 


aw -2( 2184) _ Bl (Bh) (71) 
Or|,-, OF I, (Bh) = (Bh) ) 


The rate of heat transfer can be calculated mathematically as 


(69) 


Q = re’ 


Sh=- 


00 ei Bre” Brp, gma P 
Nu = — 2 (Bh 2 L (Ph ne 
a lel ay atc |) Ca) 
5 Results 


Numerical computations were performed using the data values presented below: 
R, = 04,6 = 0.4,a = 0.35, x = 0.3,0 < Sc <10,2 < Rd, <10,0.5 < Rd, <3,2.1< Pr <21, 
1< Br <10,1.5< M <5,2<@<10,0<k <0.1,5 < Gr < 30,5 < Gc < 5. 


The results are presented in the following ways; Figs (1)-(3) represent the lipid concentration, Figs(4)-(6) 
representing the rate of mass transfer, Figs (7)-(12) representing the blood temperature profile, Figs(13)-(18) 
for rate of heat transfer, Figs(19)-(26) representing the blood velocity profile for different pertinent 
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parameter values and finally, Figs (27)-(33) showing the blood volumetric flow rate for different pertinent 
parameters. 


Fig 1 Effect on Schmidt number on Lipid Concentration with other 
values Rd, = 0.5, @ = 0.3, x = 0.3 


Fig 2 Effect on Chemical Reaction on Lipid Concentration with other values 
Sc = 0.2,@ = 0.3, x = 0.3 


Fig 3 Effect of Oscillatory Frequency on Lipid 
Concentration with other values Sc = 0.2, Rd, = 0.5, x = 0.3 
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Sh(x,t) 


Effect of Schmidt Number on Rate of Mass Transfer 
with other values Rd, = 0.5,@ = 0.3 


Sh(x,t) 


Fig 5 Effect of Chemical Reaction on Rate of Mass Transfer 
with other values Sc = 0.2. œ = 0.3 


Sh(x,t} 


Fig 6 Effect of Oscillatory Frequency on Rate of Mass 


Transfer with other values Sc = 0.2, Rd, = 0.5 
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Fig 7 Effect of Brinkman Number on Blood Temperature 
other values Sc = 0.2, @ = 0.3, Pr = 2.1, Rd, = 2, Rd, = 0.3 


Fig 8 Effect of Schmidt Number on Blood Temperature 
with other values Rd, = 0.3,@ = 0.3, Pr = 2.1, Rd, =2, Br =5 


Fig 9 Effect of Prandtl Number on Blood Temperature with 
other values Sc = 0.2,@ = 0.3, Rd, = 0.3, Rd, = 2, Br =5 
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——= Rd;=1.0 


Fig 10 Effect of Radiation Absorption on Blood Temperature 
with other values Sc = 0.2,@ = 0.3, Pr = 2.1, Rd, = 0.3, Br =5 


Fig 11 Effect of Chemical Reaction on Blood Temperature 


with other values Sc = 0.2,@ = 0.3, Pr = 2.1, Rd, =2, Br =5 


Fig 12 Effect of Oscillatory Frequency on Blood Temperature 
with other values Sc = 0.2, Rd, = 0.3, Pr = 2.1, Rd, =2, Br =5 
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Fig 13 Effect of Brinkman number on Rate of Heat Transfer 
with other values Sc = 0.2,@ = 0.3, Pr = 2.1, Rd, = 2, Rd, =0.3 


Fig 14 Effect of Schmidt Number on Rate of Heat Transfer 
with other values Rd, = 0.3,@ = 0.3, Pr = 2.1, Rd, =2, Br =5 


Fig 15 Effect of Prandtl Number on Rate of Heat Transfer 
with other values Sc = 0.2,@ = 0.3, Rd, = 0.3, Rd, = 2, Br =5 
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Nu{x,t} 


Fig 16 Effect of Radiation Absorption on Rate of Heat Transfer 
with other values Sc = 0.2,@ = 0.3, Pr = 2.1, Rd, = 0.3,Br =5 


0.2 

FoPo O A a J ----- 0.4 
z p Æ y ce W NNO J o Rd3=0.6 
=... Rd3=0.8 

——— Rd3=1.0 


Fig 17 Effect of Chemical Reaction on Rate of Heat Transfer with other 
values Sc = 0.2, @ = 0.3, Pr =2.1, Rd, =2, Br =5 


Nulx,t} 


Fig 18 Effect of Oscillatory Frequency of Heat Transfer with other 


values Sc = 0.2, Rd, = 0.3, Pr = 2.1, Rd, =2, Br =5 
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wlr, t} 


Fig 19 Effect of Brinkman number on Blood Velocity with other values Ge =15, 
Gr =15, Pr =2.1, Sc = 2, Rd, = 2, Rd, =2,@ =0.3,k = 0.05,M =1.5,x =0.3 


Fig 20 Effect of Solutal Grashof Number on Blood Velocity with other values 
Br =5,Gr =15, Pr = 2.1, Sc = 2, Rd, = 2, Rd, = 2, ø = 0.3, k = 0.05, M =1.5,x = 0.3 


Fig 20 Effect of Grashof Number on Blood Velocity with other values Br = 5, 
Ge =15, Pr =2.1, Sc = 2, Rd, =2, Rd, =2,@ =0.3,k = 0.05,M =1.5,x = 0.3 


© 2022, CAJMTCS | CENTRAL ASIAN STUDIES www.centralasianstudies.org ISSN: 2660-5309 | 47 


CENTRAL ASIAN JOURNAL OF MATHEMATICAL THEORY AND COMPUTER SCIENCES Vol: 03 Issue: 07 | Jul 2022 


Fig 21 Effect of Schmidt Number on Blood Velocity with other values Br = 5, 
Ge =15, Gr =15, Pr = 2.1, Rd, =2, Rd, = 2,0 =0.3,k =0.05,M =1.5,x =0.3 


Fig 22 Effect of Prandtl Number on Blood Velocity with other value Br = 5,Gc =15, 
Gr =15, Pr =2.1, Sc = 2, Rd, =2,Rd, =2,@ =0.3,k = 0.05,M =1.5,x =0.3 


Fig 23 Effect of Radiation Absorption on Blood Velocity with other values Br = 5, 
Ge =15, Gr =15, Pr =2.1, Sc =2, Rd, =2,@ =0.3,k =0.05,M =1.5,x =0.3 
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Fig 24 Effect of Chemical reaction on Blood Velocity with other value Br = 5, 
Ge =15, Gr =15, Pr =2.1,Se =2, Rd, =2,@ = 0.3,k =0.05,M =1.5,x =0.3 


Fig 25 Effect of Porosity on Blood Velocity with other values Br = 5,Gc = 15, 
Gr =15, Pr =2.1,Sc =2, Rd, =2, Rd, =2,@ =0.3,M =1.5,x =0.3 


Fig 26 Effect of Magnetic Field on Blood Velocity with other values 


Br =5,Ge =15,Gr =15, Pr =2.1, Sc = 2, Rd, = 2, Rd, =2,@ =0.3,k =0.05,x =0. 
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Fig 27 Effect on Brinkman number on Blood Flow Rate with other values 
Ge =15, Gr =15, Pr = 2.1, Sc =2, Rd, =2, Rd, =2,@ =0.3,k =0.05,M =1.5 


Fig 28 Effect of solutal Grashof number on Blood Flow Rate with other values 
Br =5,Gr =15, Pr =2.1,Sc = 2, Rd, =2, Rd, = 2,@ =0.3,k = 0.05,M =1.5 


Fig 29 Effect of Grashof number on Blood Flow Rate with other values 
Br =5,Ge =15, Pr =2.1, Sc =2, Rd, =2, Rd, =2,@ =0.3,k =0.05,M =1.5 
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Fig 30 Effect of Radiation Absorption on Blood Flow Rate with other values 
Br = 5,Gc =15, Gr =15, Pr =2.1, Sc =2, Rd, =2,@ =0.3,k =0.05, M =1.5 


O{x,t} 


Fig 31 Effect of chemical reaction on Blood Flow Rate with other values 
Br =5,Ge =15, Gr =15, Pr = 2.1, Sc = 2, Rd, = 2,@ = 0.3,k =0.05,M =1.5 


Fig 32 Effect of Porosity on Blood Flow Rate with other values 
Br = 5,Gc =15, Gr =15, Pr =2.1,Se = 2, Rd, =2, Rd, =2,@ =0.3,M =1.5 
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Fig 33 Effect of Magnetic Field on Blood Flow Rate with other values 


Br =5, Ge =15, Gr =15, Pr = 2.1, Sc = 2, Rd, =2, Rd, =2,@ =0.3,k =0.05 


6 Discussion 


Fig 1 shows the effect of Schmidt number on lipid concentration in a sclerotic microchannel at a specific 
location x =0.3at the timet=5. The figure indicated that the lipid concentration decreases with an increase 
in Schmidt number, this increase in Schmidt number due to low kinematic viscosity and poor diffusivity. 


Chemical reaction as result of other contributing parameter values is seen in Fig 2; the figure depicts that the 
lipid concentration decreases with an increase in chemical reaction at s specific location in the sclerotic 
artery. 


Fig 3 illustrates the effect of the oscillatory frequency on lipid concentration in the blood. This result shows 
that the lipid concentration decreases with an increasing value of the oscillatory frequency at a specific 
location x =0.3. The concentration is maximum at the centre for different oscillatory frequency and 
converged at r = 0.893 . 


The rate of mass transfer was investigated as depicted in Fig 4; the figure shows that mass transfer rate 
increases with different values of Schmidt number at different peak along the arterial channel while the other 
parameter values are Rd, = 0.5,@ = 0.3 


The effect of chemical reaction on the rate of mass transfer was investigated as seen in Fig 5, the result 
shows that the rate of mass transfer increases with an increase in chemical reaction, with other parameter 
values are Sc = 0.2, œ = 0.3 along the arterial channel. 


Fig 6 illustrates an effect of oscillatory frequency parameter change on rate of mass transfer along the 
sclerotic microchannel. This result is of the view that the mass transfer increases with an increasing value of 
the oscillatory frequency, with other values Sc = 0.2, Rd, = 0.5 along the sclerotic arterial channel. 


Fig 7 shows the effect of Brinkman number on blood temperature in a sclerotic artery at a location x = 0.3; 
the figure shows that the fluid temperature increases with an increase in Brinkman number. 


The Schmidt number effect on blood temperature was depicted in Fig 8, and this result indicates that the 
blood temperature decreases with an increase in Schmidt number in the sclerotic artery at a specific location 
x=0.3. This result is an indication that the kinematic viscosity helps in reducing the temperature rise 


Figs (9)-(12) show the effect of the Prandtl number, radiation absorption, chemical reaction parameter and 
oscillatory frequency parameter on the blood temperature. These results indicate that the blood temperature 
decreases with an increase in Prandtl number, radiation absorption, chemical reaction and oscillatory 
frequency at a specific location x = 0.3in a sclerotic artery with treatment R, = 4 in timer =5. 
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Figs (13)-(18) illustrate the effects of Brinkman number, Schmidt number, Prandtl number, radiation 
absorption, chemical reaction and oscillatory frequency on the rate of heat transfer along sclerotic arterial 
channel. The results indicate that the rate of heat transfer increase with different values of the 
aforementioned parameters along the vessel at the treatment R, = 0.4 , growth rate a = 0.34 . 


The Brinkman number effect on blood flow was investigated as depicted by Fig 19; the figure shows that the 
velocity increases for an increase in Brinkman number with other parameter values 
Gr =15, Pr =2.1,Sc =2, Rd, =2, Rd, =2,@ =0.3,k =0.05,M =1.5 at a specific location in the sclerotic 


microchannel. 


Thermal Grashof number effect on blood velocity was investigated and result shown in Fig 20; the result 
shows that blood velocity increases with an increase in Grashof number at a specific location x =0.3, with 
other parameter values such as Gc=15,Pr=2.1,Sc=2, Rd, =2, Rd, =2,@=0.3,k =0.05,M =1.5, 


contributing to this flow behavior. 


Fig 21 depicts the effect of Schmidt number on blood velocity in a sclerotic microchannel at a specific 
location x = 0.3. The figure shows that the blood velocity decreases with an increase in Schmidt number and 
the velocity decrease is due to low diffusivity. 


The effect on Prandtl number on blood velocity was investigated as seen in Fig 22; the result shows that the 
velocity decreases for an increasing value of the Prandtl number. 

Fig 23 illustrates the effect of radiation absorption on blood velocity in a sclerotic microchannel artery. It is 
seen that the blood velocity reduces as radiation absorption increases from Rd, = 2,4,6,8,10at x =0.3. 


The chemical reaction as a result of the treatment and other contributing factors denoted with the parameters 
Ge =15,Gr =15, Pr =2.1,Sc =2, Rd, =2,@ =0.3,k =0.05, M =1.5,tends to reduce the velocity of the fluid 


in a sclerotic vessel as seen in Fig 24. The figure shows that blood velocity decreases for an increasing value 
of chemical reaction. 


The study investigated the blood velocity at different of porosity in a sclerotic vessel as depicted by Fig 25. 
The figure shows that the blood velocity increases at different level of porosity, and it’s maximum at the 
centre line of the artery for different porosity as well. 


Fig 26 illustrates the effect of magnetic field on blood velocity at a specific location of x=0.3in the 
sclerotic microchannel att =5. In this figure, the axial velocity reduces with an increasing magnetic field 
parameter. This occurs due to the interaction of magnetic field with blood which produces body force called 
the Lorentz force, which impede the motion of blood in the vessel. The maximum velocity was observed at 
the centre line of the artery as depicted in the figure. 


The research focuses on the effect of Brinkman number on blood flow; however, we investigated the effect 
of the Brinkman number on volumetric flow rate as seen in Fig 27. This figure indicates that the volumetric 
flow rate increases for an increase in Brinkman number. This useful flow rate behavior is made possible with 
the contribution of other pertinent parameters such asGc=15,Gr=15,Pr=2.1,Sc = 2, Rd, =2, 


Rd, =2,@=0.3, k =0.05, M =1.5. 
Fig 28 indicates that the flow rate increases for different values of the solutal Grashof. The result is of the 
view that for us to increase the volumetric flow of viscous fluid as blood, we need to increase the solutal 


Grashof number which is due to the lipid concentration i the fluid. However, the flow is not complete 
without the contribution of other parameter values such as 


Br =5,Gr =15, Pr =2.1, Sc =2, Rd, =2, Rd, = 2,@=0.3,k =0.05, M =1.5. 


The thermal Grashof number was investigated along the arterial channel with other parameters values such 
as Br =5,Ge =15, Pr = 2.1, Sc =2, Rd, =2, Rd, =2,@=0.3, k =0.05, 
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M =1.5. This result is of the view that, the flow rate increases for different values of the thermal Grashof 
number as shown in Fig 29. 


The radiation absorption was investigated and the result displayed in Fig 30, which shows that the flow rate 
decreases for an increase in radiation absorption with other parameters values 
Br = 5,Gc =15,Gr =15, Pr = 2.1, Sc = 2, Rd, =2,@=0.3,k =0.05,M =1.5 contributing to that flow behavior. 


Fig 31 illustrates the blood volumetric flow rate for different values of the chemical reaction. This figure is 
of the views that increase in chemical reaction triggers a decrease in volumetric flow rate. 


Fig 32 shows the blood flow rate for different values of the porosity in the blood vessel at t = 5 and the other 
contributing values are Br = 5,Gc =15, Gr =15, Pr =2.1,Sc =2, Rd, =2, Rd, =2,@=0.3,M =1.5. Magnetic 
field effect on flow rate was investigated as seen in Fig 33; the figure shows that the volumetric flow rate of 
blood decreases with increasing values of magnetic field along the arterial channel with treatment and 
growth rate of R, =0.4 and a = 0.34 respectively. 


7 Conclusion 


The radiation-brinkman number effects on blood flow in a porous atherosclerotic microchannel with the 
presence of a magnetic field, growth rate, and treatment. Using the Laplace method, the nonlinear partial 
differential equations with the effective governing parameters were examined. Our results were revealed as 
follows: 


a. The axial velocity and volumetric flow decrease with the increase in Schmidt number, Prandtl number, 
radiation absorption, chemical reaction, magnetic field parameter. However, as the Brinkman number, 
solutal Grashof number, thermal Grashof number, and porosity parameter increase, so do the flow 
characteristics. 


b. The cholesterol concentration decreases with an increase in Schmidt number, chemical reaction and 
oscillatory frequency, while the mass transfer rate increases with an increase in Schmidt number, 
chemical reaction and oscillatory frequency. 


c. The blood temperature rises with an increase in Brinkman number but decreases with an increase in 
Schmidt number, Prandtl number, radiation absorption, chemical reaction, and oscillatory frequency, 
while the rate of heat transfers increases with an increase in Brinkman number, Schmidt number, Prandtl 
number, radiation absorption, chemical reaction, and oscillatory frequency. 


Nomenclature 
x Dimensional coordinate along the channel 
r Dimensional coordinate perpendicular to the channel 


R Radius of an abnormal channel 


Ro Radius of normal channel 

P Systolic pressure 

Rd, Radiation absorption parameter 
Rd, Chemical reaction parameter 
k,,, Blood thermal conductivity 


Dimensional velocity profile 


Dimensionless velocity profile 
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Wo Perturbed velocity profile 


C Dimensional lipid particle concentration 
C, Far field cholesterol particle concentration 
Cop The specific heat capacity of blood 

f Dimensionless time 


T Temperature of the fluid 


T. Far field temperature 
Temperature at the wall 
B, Magnetic intensity 


M Magnetic field parameter 


a Growth rate of LDL-cholesterol 
Greek Symbols 
V, Kinematic viscosity of blood 


L, Dynamic viscosity of blood 


Pr Prandtl number for blood 


Acceleration due to gravity 


Dimensional height of stenosis 


D, Electrical conductivity 
A” Length of stenosis 
a) Oscillatory frequency 


Dimensionless blood temperature 


0 
o Dimensionless cholesterol particle concentration 
0 Dimensionless wall temperature 

0, Perturbed blood temperature profile 

?, Dimensionless wall lipid concentration 

ho Perturbed lipid concentration 


P, Blood density 


Subscripts 

w Wall 

b Blood 

e Electrical 
T Thermal 
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Far field 


MMDARG Mathematical Modeling and Data Analytic Research Group 
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